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[i]  During  the  main  phase  of  the  6  April  2000  magnetic  storm  with  Dst  =  —300  nT,  four 
Defense  Meteorological  Satellite  Program  (DMSP)  satellites  encountered  intense  sheets  of 
field-aligned  currents  (FACs).  Their  magnetic  perturbations  were  >1300  nT, 
corresponding  to  integrated  currents  |J|||  >  1  A/m.  The  FACs  appeared  in  both  the  evening 
and  dawn  magnetic  local  time  sectors.  They  had  relatively  fast  rise  times  (~5  min),  lasted 
for  ~0.5  hours,  and  were  associated  with  widespread  reconfigurations  of  plasma  in  the 
near-Earth  magnetotail.  The  largest  magnetic  and  related  electric  field  perturbations 
occurred  at  magnetic  latitudes  <60°.  Magnetometer  measurements  from  DMSP  satellites 
show  repeated  episodes  of  similarly  large  FACs  late  in  the  main  phase  of  this  and 
other  superstorms.  Poynting  flux  calculations  indicate  that  a  few  percent  of  the  total  ring 
current  energy  is  dissipated  as  Joule  heat  in  the  midlatitude  ionosphere  during  each  of 
these  events.  A  survey  of  ground  magnetometers  at  auroral  and  middle  latitudes  found 
perturbations  typically  <200  nT,  incommensurate  with  magnetic  observations  at  the 
altitude  of  the  DMSP  satellites.  The  small  ground  responses  reflect  weak  ionospheric  Hall 
currents  and  provide  an  empirical  validation  of  the  theorem  derived  by  Fukushima  [1976]. 

Height-integrated  Pedersen  conductances  (EP)  calculated  with  observed  precipitating 
electron  fluxes  have  significantly  lower  values  than  those  estimated  from  Ohm’s  law  using 
simultaneously  measured  electric  and  magnetic  field  variations.  This  discrepancy  suggests 
that  under  some  highly  stressed  conditions  the  neglect  of  contributions  from  low  , 
energy  electrons  and  precipitating  ions  can  lead  to  serious  underestimates  of  EP  and  a 
consequent  misunderstanding  of  the  magnetosphere-ionosphere  circuit.  Index  Terms: 

2409  Ionosphere:  Current  systems  (2708);  2431  Ionosphere:  Ionosphere/magnetosphere  interactions  (2736); 

2788  Magnetospheric  Physics:  Storms  and  substorms;  2716  Magnetospheric  Physics:  Energetic  particles, 
precipitating;  KEYWORDS:  field-aligned  currents,  superstorms,  measured  conductances,  DMSP 
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1.  Introduction 

[2]  We  report  the  first  observations  of  intense  (~1  A/m) 
but  transient  sheets  of  field-aligned  currents  (FACs)  that 
couple  the  central  plasma  sheet  with  the  auroral  ionosphere 
near  the  maximum  epochs  of  magnetic  superstorms.  Max¬ 
imum  magnetic  perturbations  associated  with  the  FACs 
normally  develop  at  middle  (55° -60°)  magnetic  latitudes. 
Our  attention  was  drawn  to  these  phenomena  while  com¬ 
paring  ground  magnetometer  responses  to  measurements 
by  particle  and  field  sensors  on  Defense  Meteorological 
Satellite  Program  (DMSP)  satellites  during  a  4-hour  period 
on  6  April  2000.  Six  empirical  aspects  of  the  data  appeared 
to  be  unusual: 
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[3]  1 .  Although  four  DMSP  satellites  traversed  the  FACs 
at  about  the  same  time  in  the  evening  and  dawn  magnetic 
local  time  (MLT)  sectors,  no  commensurate  magnetic  per¬ 
turbations  were  detected  at  ground  stations  beneath  or 
conjugate  to  them.  This  null  result  was  confirmed  in  a 
survey  of  data  from  many  magnetometer  stations  at  auroral 
and  middle  latitudes. 

[4]  2.  Pedersen  conductances  (EP)  calculated  using 
Ohm’s  law  differed  from  estimates  based  on  energy  fluxes 
of  precipitating  electrons. 

[5]  3.  In  the  evening  sector,  precipitating  ion  fluxes 
dominate  over  precipitating  electrons  during  the  interval 
of  downward  FAC. 

[6]  4.  As  the  evening-sector  current  system  decayed,  a 
new  upward  FAC  developed  near  the  equatorward  boundary 
of  auroral  electron  precipitation.  Ion  fluxes  decreased 
dramatically,  suggesting  that  upward  drifting,  low-energy 
electrons  then  carried  the  ~1  A/m  downward  current. 
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[7]  5.  A  significant  fraction  of  the  total  ring  current 
energy  was  dissipated  as  Joule  heat  in  the  midlatitude 
ionosphere. 

[8]  6.  Four  distinct  episodes  of  intense  FACs  occurred 
during  the  storm’s  maximum  epoch,  each  following  a 
pseudo-breakup  event,  apparent  in  traces  of  the  auroral 
electrojet  (AE)  index. 

[9]  Both  individually  and  collectively,  these  unusual 
DMSP  observations  suggest  that  under  highly  stressed 
conditions,  magnetosphere-ionosphere  (M-I)  coupling  pro¬ 
ceeds  in  ways  not  fully  grasped  in  current  models.  The 
purpose  of  this  paper  is  twofold.  First,  we  describe  relevant 
DMSP  measurements  in  the  geophysical  context  of  a  severe 
magnetic  storm.  Second,  we  relate  DMSP  data  to  the  current 
understanding  of  M-I  coupling  derived  from  theoretical  and 
observational  considerations.  We  thus  identify  gaps  between 
DMSP  observations  and  existing  models  as  indicators  of 
new  paths  for  future  exploration. 

[10]  From  today’s  perspective  it  is  difficult  to  realize  that 
the  existence  of  FACs  was  a  hotly  contested  issue  just 
30  years  ago  [Egeland,  1984;  Dessler,  1984].  The  existence 
debate  effectively  ended  with  publications  by  Iijima  and 
Potemra  [1976,  1978],  showing  that  global  distributions  of 
FACs  exist  and  follow  relatively  simple  patterns  to  which 
they  assigned  the  durable  monikers  of  Region  1  (Rl)  and 
Region  2  (R2).  Rl  currents  are  located  near  the  poleward 
boundary  of  the  auroral  oval,  flowing  into  the  ionosphere  on 
the  dawn  side  of  the  auroral  oval  and  out  on  the  dusk  side. 
R2  currents  are  equatorward  of  Rl  and  have  opposite 
polarities,  flowing  into  the  ionosphere  on  the  duskside 
and  out  on  the  dawnside. 

[11]  Two  key  papers  of  the  FAC  debate  are  particularly 
relevant  to  the  present  study.  The  first  concerns  physical 
requirements/implications  for  FACs  as  agents  responsible 
for  M-I  coupling.  Vasyliunas  [1970]  used  force  balance  and 
current  continuity  to  show  that  FACs  are  consequences  of 
misalignments  between  plasma  pressure  gradients  and  flux- 
tube  volume  gradients  in  the  magnetosphere.  Ohm’s  law 
dictates  that  FACs  control  distributions  of  electric  potential 
in  the  ionosphere.  Indirectly,  through  magnetic  mapping, 
they  also  control  electric  field  distributions  in  the  magneto¬ 
sphere.  The  second  paper  concerns  ground  signatures  of 
FACs.  Fukushima  [1976]  developed  a  theorem  demonstrat¬ 
ing  that  if  FACs  close  through  poloidal  currents  flowing  in  a 
uniformly  conducting  ionosphere,  magnetic  perturbations  at 
ground  level  exactly  cancel  and  produce  no  ground  signa¬ 
tures.  We  refer  to  this  result  as  Fukushima’s  theorem. 

[12]  The  Rice  Convection  Model  (RCM)  implements 
Vasyliunas’  prescription  [Harel  et  al.,  1981].  RCM  only 
simulates  the  dynamics  of  the  inner  magnetosphere,  treating 
Rl  currents  and  coupling  to  interplanetary  space  as  outside 
its  boundary  of  concern.  It  assumes  that  plasma  compresses 
adiabatically  while  calculating  the  growth  of  the  ring  current, 
the  formation  of  R2  currents,  and  consequent  redistributions 
of  electric  fields  in  the  ionosphere/inner  magnetosphere. 
However,  Erickson  and  Wolf[  1980]  pointed  out  that  adia¬ 
batic  compression  leads  to  insupportably  high  plasma  pres¬ 
sures.  Kivelson  and  Spence  [1988]  suggested  that  ion  loss 
via  gradient-curvature  drift  to  the  magnetopause  lowers 
plasma  pressure.  Borovsky  et  al.  [1998]  empirically  demon¬ 
strated  that  significant  departure  from  adiabatic  transport 
occurs  at  radial  distances  between  15  and  6.6  RE  in  the 


magnetotail.  They  compared  several  possible  mechanisms, 
concluding  that  auroral  ion  precipitation  is  responsible  for  a 
large  fraction  of  the  mass  loss.  We  note  that  a  recent  RCM 
simulation  of  the  June  1991  magnetic  storm  [Gamer  et  al., 
2004]  reproduced  electric  field  distributions  observed  by  the 
Combined  Release  and  Radiation  Effects  Satellite  (CRRES) 
in  the  inner  magnetosphere  and  DMSP  satellites  in  the 
ionosphere  [Burke etal.,  1998].  However,  RCM  consistently 
predicted  more  plasma  pressure  than  was  measured  in  the 
ring  current  population  [Garner  et  al.,  2004],  While  de  facto, 
pressure  balance  is  maintained  in  the  magnetosphere,  details 
concerning  how  this  is  achieved  remain  unclear. 

[13]  Richmond  and  Kamide  [1988]  developed  an  empir¬ 
ical  approach  to  specify  FAC  and  ionospheric  potential 
distributions  that  grew  into  detailed  procedures  for  integrat¬ 
ing  ground  and  space  data  called  assimilative  mapping  of 
ionospheric  electrodynamics  (AMIE)  [Richmond,  1992]. 
Magnetic  perturbations  observed  by  ground  stations  at  high 
and  middle  latitudes  are  inverted  to  specify  equivalent 
ionospheric  currents.  Raeder  et  al.  [2001]  argued  that 
translating  equivalent  currents  into  realistic  electric  field 
estimates  requires  precise  knowledge  about  the  distribution 
of  ionospheric  conductances.  To  varying  degrees  of  success, 
conductances  can  be  derived  from  knowledge  of  solar 
ultraviolet  radiances  [Wallis  and  Budzinski,  1981]  and  the 
fluxes  of  energetic  particles  precipitating  from  the  magne¬ 
tosphere  [Robinson  et  al.,  1987].  However,  Kamide  et  al. 
[1981]  warned  that  perturbations  produced  by  FACs  and 
poloidal  closing  currents  exactly  cancel  on  the  ground. 
AMIE  is  thus  vulnerable  to  underestimating  FAC  intensities 
and  ionospheric  electric  fields. 

[14]  M-I  coupling  involves  exchanges  of  charged  par¬ 
ticles  as  well  as  electromagnetic  energy.  The  morphologies 
and  origins  of  energetic  particle  precipitation  into  the  polar 
cap  and  nightside  auroral  oval  are  relevant  to  this  report. 
Winningham  and  Heikkila  [1974]  showed  that  the  polar  cap 
is  characterized  by  nearly  uniform  precipitation  of  electrons 
with  E  <  100  eV,  called  polar  rain.  Winningham  et  al. 
[1975]  distinguished  between  two  types  of  auroral  electron 
precipitation  referred  to  as  boundary  plasma  sheet  (BPS) 
and  central  plasma  sheet  (CPS)  precipitation.  BPS  precip¬ 
itation  occurs  in  the  poleward  part  of  the  nightside  oval  and 
is  highly  structured  in  energy-versus-time  spectrograms. 
CPS  precipitation  occurs  in  the  equatorward  part  of  the 
oval  and  varies  slowly  ,  with  latitude.  DeCoster  and  Frank 
[1979]  used  a  similar  nomenclature  to  describe  magneto- 
spheric  populations  observed  in  the  “plasma  sheet  bound¬ 
ary  layer”  (PSBL)  and  “central  plasma  sheet,”  leading  to 
confusion  about  relationships  between  BPS  and  PSBL. 
Galperin  and  Feldstein  [1991]  showed  that  the  PSBL 
covered  a  much  narrower  range  of  magnetic  latitudes  than 
the  BPS.  Episodically,  DMSP  observes  PSBL  signatures 
near  the  poleward  boundary  of  the  auroral  oval  as  ~  1 0  wide 
regions  of  intense  electric  fields  and  velocity-dispersed  ion 
fluxes  [Burke  et  al.,  1994], 

[15]  Recently,  Newell  et  al.  [1996]  sorted  ion  and  electron 
fluxes  by  latitude  in  the  auroral  oval  and  related  the  bound¬ 
aries  to  transitions  in  the  conjugate  magnetosphere.  This 
extends  work  by  Gussenhoven  et  al.  [1983]  who  interpreted 
the  equatorward  boundary  of  auroral  electron  precipitation 
as  mapping  to  the  zero-energy  Alfven  boundary  in  the 
magnetosphere.  Four  of  the  boundaries  designated  by 
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Newell  and  coworkers  are  useful  here.  Proceeding  from  low 
to  high  latitude,  these  are  (1)  the  equatorward  boundary  of 
auroral  precipitation  Bl,  equivalent  to  the  instantaneous 
zero-energy  Alfven  boundary,  (2)  entry  into  the  main 
plasma  sheet  B2e,  where  average  electron  energies  cease 
to  increase  with  magnetic  latitude,  (3)  the  ion  isotropy 
boundary  B2i,  marking  the  transition  between  quasi-dipolar 
and  stretched  field  lines,  and  (4)  the  poleward  boundary  of 
the  main  oval  B5. 

[i6]  Burke  et  al.  [1998]  compared  distributions  of  ener¬ 
getic  particles  and  electric  fields  observed  by  CRRES  in  the 
inner  magnetosphere  with  those  detected  by  DMSP  F8  in 
the  ionosphere  during  the  magnetic  superstorm  of  4-6  June 
1991.  In  this  study  the  phrase  “penetration  electric  fields” 
referred  to  electric  fields  in  the  magnetosphere  earthward  of 
the  electron  plasma  sheet  boundary,  which  is  magnetically 
conjugate  to  Bl.  As  the  storm’s  main  phase  progressed,  Bl 
moved  equatorward  and  stabilized  at  ~53°  magnetic  lati¬ 
tude  (MLat).  In  a  dipole  mapping  this  corresponds  to  L  «  3. 
However,  Figure  6  of  Burke  et  al.  [1998]  shows  that 
throughout  the  main  phase  CRRES  crossed  the  inner  edge 
of  the  electron  plasma  sheet  beyond  L  =  6.  The  difference 
between  the  dipole  and  real  mapping  of  Bl  is  caused  by  the 
inflationary  effects  of  the  ring  current  on  magnetic  fields 
near  the  magnetosphere’s  equatorial  plane  [Tsyganenko  and 
Stern,  1996]. 

[n]  The  following  section  describes  orbital  character¬ 
istics  and  scientific  payloads  of  DMSP  satellites  and  the 
spacecraft-centered  coordinate  system  used  for  data  pre¬ 
sentation.  Before  presenting  DMSP  observations,  we  sum¬ 
marize  our  techniques  for  extracting  M-I  coupling 
parameters  from  combined  particle  and  field  data.  The 
observations  section  includes  an  overview  of  the  6  April 
2000  “superstorm,”  a  term  used  to  describe  magnetic 
disturbances  with  Dst  <  -200  nT.  A  geophysical  context 
for  understanding  DMSP  observations  is  established 
through  consideration  of  solar  wind  and  interplanetary 
magnetic  field  (IMF)  dynamics  and  geospace  responses 
represented  by  variations  of  selected  geomagnetic  indices 
and  auroral  boundaries  [Newell  et  al.,  1996],  The  rest  of 
the  observations  section  presents  the  DMSP  data  required 
to  support  the  six  empirical  results  enumerated  in  the 
introduction.  The  discussion  section  addresses  (1)  some 
implications  of  DMSP  measurements  for  understanding  the 
dynamics  of  plasma  pressure  distributions  in  the  storm¬ 
time  magnetosphere,  (2)  applications  of  Ohm’s  law  in  the 
ionosphere,  and  (3)  validation  of  Fukushima’s  theorem. 

2.  Instrumentation 

[is]  DMSP  satellites  are  three-axis  stabilized  spacecraft 
that  fly  in  circular,  sun-synchronous,  polar  (inclination 
98.7°)  orbits  at  an  altitude  of  ~840  km.  The  geographic 
local  times  of  the  orbits  are  either  near  the  1800-0600  (FI  1, 
F13)  or  2100-0900  (F12,  14,  15)  local  time  meridians. 
Offsets  between  the  geographic  and  geomagnetic  poles 
allow  DMSP  satellites  to  sample  wide  ranges  of  magnetic 
local  times  over  the  course  of  a  day.  The  ascending  nodes  of 
DMSP  orbits  are  on  the  dusk  side  of  the  Earth.  Thus  the 
satellites  move  toward  the  northwest  (southeast)  in  the 
evening  (morning)  sector.  Each  satellite  carries  a  suite  of 
sensors  to  measure  (1)  fluxes  of  auroral  particles  (Special 


Sensor  J4,  SSJ4)  and  (2)  densities,  temperatures,  and  drift 
motions  of  ionospheric  ions  and  electrons  (Special  Sensor 
for  Ions  Electrons  and  Scintillations,  SSIES).  All  of  the 
satellites  except  Fll  carry  magnetometers  (Special  Sensor 
for  Magnetic  Fields,  SSM)  to  monitor  perturbations  of  the 
Earth’s  magnetic  field. 

[19]  Identical  SSJ4  sensors  are  mounted  on  the  top  surfaces 
of  DMSP  satellites  to  measure  fluxes  of  down-coming 
electrons  and  ions  in  19  logarithmically  spaced  energy  steps 
between  approximately  30  eV  and  30  keV  [Hardy  et  al., 
1984],  A  ftill  spectrum  is  compiled  every  second.  Unfortu¬ 
nately,  fluxes  of  ions  with  energies  below  1  keV  are  unavail¬ 
able  from  the  FI 5  satellite.  SSIES  consists  of  spherical 
Langmuir  probes  mounted  on  0.8  m  booms  to  measure 
densities  and  temperatures  of  ambient  electrons,  and  three 
different  sensors  mounted  on  a  conducting  plate  facing  the 
ram  direction.  These  are  (1)  ion  traps  to  measure  the  total  ion 
densities,  (2)  ion  drift  meters  to  measure  horizontal  (VH)  and 
vertical  (Vv)  cross-track  components  of  plasma  drifts,  and  (3) 
retarding  potential  analyzers  (RPA)  to  measure  ion  temper¬ 
atures  and  in-track  components  of  plasma  drift  V||  [Rich  and 
Hairston,  1994],  Shapes  of  RPA  current- voltage  sweeps  can 
be  used  to  determine  the  masses  of  contributing  ions 
species.  Because  of  relatively  small  uncertainties  about 
spacecraft  potentials,  RPA  estimates  of  Vpar  often  have  large 
error  bars  and  are  not  used  in  this  study.  Although  the  drift 
meters  measure  VH  six  times  per  second,  data  are  presented 
here  as  4-s  averaged  values.  Problems  with  a  noisy  tape 
recorder  limit  the  usefulness  of  data  from  Fll. 

[20]  SSM  sensors  are  triaxial  fluxgate  magnetometers  that 
are  mounted  on  the  bodies  of  the  F12-F14  spacecraft.  The 
magnetometer  on  FI 5  is  mounted  on  a  5-m  boom,  reducing 
its  susceptibility  to  spacecraft-generated  contamination. 
Magnetic  vectors  are  sampled  12  times  per  second.  One- 
second  average  values  are  calculated  as  8B  =  Bmeas  -  BIGRF, 
the  differences  between  measured  and  International  Geo¬ 
magnetic  Reference  Field  (IGRF)  values  of  magnetic  fields 
at  the  spacecraft  locations.  Data  are  presented  as  8BX,  8BY, 
6BZ  in  spacecraft-centered  coordinates,  as  illustrated  in 
Figure  1 .  The  X-  and  Y-axes  point  toward  spacecraft  nadir 
and  along  the  velocity  vector,  respectively.  The  Z-axis 
completes  the  right-hand  system.  Figure  1  shows  projec¬ 
tions  onto  an  invariant  latitude-versus-MLT  grid  of  four 
DMSP  trajectories  near  2020  UT  on  6  April  2000.  FI 3  and 
FI 4  (FI 2  and  FI 5)  were  moving  equatorward  in  the 
morning  (evening)  sector  above  the  northern  (southern) 
ionosphere.  The  Z-axis  generally  points  in  the  antisunward 
direction,  as  shown  in  Figure  1.  Nearly  simultaneous 
observations  from  the  SSM,  SSIES,  and  SSJ4  sensors  on 
the  various  DMSP  spacecraft  help  corroborate  and  interpret 
electrodynamic  measurements  in  the  storm’s  main  phase. 

[21]  To  assist  the  interpretation  of  DMSP  measurements 
we  also  examined  other  data  sources.  These  include: 
(1)  interplanetary  plasma  and  magnetic  field  measurements 
from  the  Wind  and  Advanced  Composition  Explorer  (ACE) 
satellites,  (2)  magnetic  fields  measured  at  a  large  number  of 
ground  stations  distributed  around  the  Earth,  and  (3) 
observations  of  energetic  (keV  and  MeV)  particle  fluxes 
by  four  Los  Alamos  National  Laboratory  (LANL)  satellites 
and  magnetic  field  by  two  Geostationary  Operational 
Environmental  Satellite  (GOES)  satellites  at  geostationary 
altitude.  The  LANL  satellites  detected  fluctuations  in  keV 


3  of  15 


A06303 


HUANG  AND  BURKE:  FIELD-ALIGNED  CURRENTS  DURING  A  SUPERSTORM 


A06303 


06:00 


18:00 

Figure  1.  Spacecraft-centered  coordinate  system  used  to 
analyze  SSM  data.  Trajectories  of  four  DMSP  satellites  are 
projected  onto  an  invariant  latitude  versus  magnetic  local 
time  grid.  X  and  Y  are  directed  toward  nadir  and  along  the 
spacecraft  velocity  vector,  respectively.  Z  generally  points 
toward  the  antisunward  direction.  Black  circles  mark 
locations  on  the  trajectories  where  maximum  perturbations 
in  6BZ  were  detected. 

particle  fluxes  that  intensified  with  increased  geomagnetic 
activity.  The  observed  variations  showed  no  obvious  corre¬ 
lations  with  measurements  from  DMSP  or  each  other.  The 
GOES  satellites  were  both  on  the  dayside  during  the  period 
of  interest  and  provide  no  information  about  nightside 
electrodynamics. 

3.  Data  Analysis  Techniques 

[22]  This  paper  integrates  nearly  simultaneous  measure¬ 
ments  from  three  sensors  on  four  DMSP  satellites  to  draw 
inferences  about  M-I  coupling  under  highly  stressed  con¬ 
ditions.  Techniques  used  to  draw  geophysical  inferences 
from  the  measurements  by  each  sensor  are  distributed  in  the 
literature  and  are  very  familiar  to  specialized  analysts. 
However,  it  appears  useful  for  the  benefit  of  more  general 
readers  to  summarize  here  the  analysis  techniques  and  their 
interrelationships.  This  section  treats  characteristics  of  data 
streams,  in  spacecraft-centered  coordinates,  expected  when 
DMSP  satellites  cross  a  pair  of  quasi-infinite  FAC  sheets  at 
normal  incidence.  Methods  for  identifying  breakdowns  of 
the  normal  incidence  and/or  infinite  current-sheet  approx¬ 
imations  are  briefly  indicated. 

[23]  The  two  governing  equations  for  analyzing  com¬ 
bined  SSM  and  SSIES  data  are  Ampere’s  law  and  current 
continuity: 


j||  —  O/FoM^  x  »]| 

(1) 

3j||/<?s  =  -Vx  Mx, 

(2) 

where  s  represents  an  infinitesimal  distance  along  a 
magnetic  field  line.  In  the  ionosphere  the  current  density 
perpendicular  to  the  Earth’s  magnetic  field  BE  is 

jx  =  [dPE  -  cth(E  x  b)],  (3) 


where  o>  and  crH  are  the  Pedersen  and  Hall  conductivity, 
respectively;  b  is  a  unit  vector  along  Be-  Integrating 
equation  (2)  from  the  satellite  location  to  the  bottom  of  the 
ionosphere  gives 

j|l  =  Vx.Ii=V1.[£PE  —  EH(Exb)],  (4) 

where  XP  and  £H  are  the  Pedersen  and  Hall  conductances, 
respectively. 

[24]  Applied  to  a  DMSP  satellite  crossing  an  infinite 
field-aligned  current  sheet  at  normal  incidence  Vx  — >  dy, 
and  equation  (1)  reduces  to 

j|,  =  (1/Mo)  [SffiBz].  (1') 

By  definition,  FACs  into  the  ionosphere  flow  in  the  +X 
direction.  Equation  (l')  indicates  that  in  plots  of  6BZ  versus 
time,  FACs  into  the  ionosphere  are  characterized  by  positive 
slopes.  Since  DMSP  satellites  move  in  the  +Y  direction  at 
speeds  of  ~7.5  km/s,  a  FAC  density  of  1  pA/m2 
corresponds  to  a  slope  of  ~9.4  nT/s. 

[25]  A  useful  quantity  for  estimating  the  total  current  into 
or  out  of  the  ionosphere  is  the  integrated  (A/m)  field-aligned 
current  density  J||  =  J  j||  dY.  Integration  of  equation  (l') 
along  a  trajectory  segment  gives 

J,,  =  (1/p,)  [A6BZ],  (5) 

where  the  symbol  A  represents  the  difference  between  6BZ 
values  at  the  applied  limits  of  integration.  Jj|  =  1  A/m 
corresponds  to  a  magnetic  perturbation  of  1256  nT. 

[26]  In  our  spacecraft-centered  coordinate  system,  equa¬ 
tion  (4)  becomes 

j||  =  0yPp  Ey  -  £h  Ez].  (4') 

Note  that  for  a  normal  incidence  crossing,  Ez  is  the  electric 
field  component  tangential  to  the. sheet  direction.  Faraday’s 
law  requires  that  to  be  constant.  Combining  equations  (l') 
and  (4')  gives 

9y[8Bz  -  p,(Ep  Ey  -  Eh  Ez)]  =  0.  (6) 

In  analyzing  combined  SSM  and  SSIES  measurements,  it  is 
often  useful  to  identify  segments  in  data  streams  where 
variations  of  6BZ  and  EY  are  highly  correlated.  Equation  (6) 
indicates  that  during  such  data-acquisition  segments  the 
conductance  gradients  are  weak.  This  judgment  can  be 
tested  against  the  properties  of  particle  fluxes  measured  by 
the  onboard  SSJ4  sensor  during  the  same  interval  [ Robinson 
et  al.,  1987].  Within  the  regions  of  weak  conductance 
gradients 

EP  «  (1/p,)  [A8Bz/AEy],  (7) 

Our  analyses  of  DMSP  measurements  also  identify  particle 
populations  responsible  for  jy.  SSJ4  sensors  monitor  ion/ 
electron  fluxes  moving  in  the  +X  direction.  Differential 
fluxes  measured  by  SSJ4  may  be  integrated  over  energy 
(30  eV  to  30  keV)  to  determine  number  fluxes  4?N,s(tp,  \j /), 
in  particles  per  (cm2  s  sr),  for  each  species,  s.  Here  ip  and 
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\| /  represent  the  gyrophase  and  pitch  angles,  respectively. 
Assuming  that  both  species  are  gyrotropic, 

j|l  =  2it  q  j  [$N,iOp,  Of)  -  $M,e((p,  a)]Cos  \j/  Sin  y  d \\i,  (8) 

where  q  represents  the  elemental  unit  of  charge  1.62  x 
10~19  C.  Consider  a  sheet  of  upward  FAC  and  assume 
that  both  precipitating  particle  distributions  are  iso¬ 
tropic  over  the  down-coming  hemisphere,  then 

j  ||  ~  ir  q  [$N,i  (V,  <*)  -  $N,e  (ip,  a)] .  (9) 

Equation  (9)  is  valid  as  long  as  backscattered  fluxes  are 
much  less  than  the  incident  fluxes.  Experience  teaches 
that  in  most  geophysical  cases  upward  and  downward 
FACs  are  carried  by  electrons  moving  into  and  out  of 
the  ionosphere,  respectively.  During  the  6  April  2000 
magnetic  storm,  DMSP  detected  an  unusual  case  in 
which  an  intense  downward  FAC  was  carried  by 
precipitating  energetic  ions  (empirical  result  3). 

[27]  The  final  physical  quantity  used  as  an  analysis  tool  in 
this  paper  is  the  field-aligned  Poynting  vector  Sy  (W/m2)  = 
(l/p0)  [E  x  8B]||  «  (EY  6Bz/po)  x.  For  reference,  Ey  and 
6BZ  perturbations  of  1  mV/m  and  1  nT  correspond  to  S||  « 
0.8  pW/m2.  On  the  time  scales  of  interest,  Poynting’s 
theorem  relates  the  divergence  in  S||  to  the  Joule  heat 
dissipated  in  the  ionosphere  beneath  the  spacecraft 

d  Sy/ds  =  — jx  •  E  =  -crP  E2.  (10) 

Integrating  equation  (10)  from  the  satellite  altitude  to  the 
bottom  of  the  ionosphere  gives  ASy  «  EP  Ey.  The  symbol 
AS  ||  represents  the  difference  in  electromagnetic  energy  flux 
carried  by  incident  and  reflected  Alfven  waves.  In  the 
discussion  section  we  show  that  ASy  is  the  physical  quantity 
inferred  from  measured  values  of  EYand  6BZ.  Equation  (10) 
thus  indicates  that  the  total  Joule  heat  dissipated  in  the 
ionosphere  can  be  estimated  from  combined  Ey  and  6BZ 
measurements  with  no  prior  knowledge  of  Ep.  The  total  rate 
of  electromagnetic  energy  input  to  the  ionosphere  associated 
with  sheets  of  FAC  is 

nn(W)  «  y  ASy  dY  dZ.  (11) 

While  we  have  no  direct  knowledge  of  the  current  sheet’s 
extent  in  the  ±Z  directions,  we  can  integrate  equation  (11) 
numerically  along  spacecraft  trajectories  to  determine 
electromagnetic  power  (and  Joule  heat)  in  Watts  per  meter 
transferred  from  the  magnetosphere  to  the  ionosphere. 

4.  Observations 

[28]  Particle  and  field  measurements  by  sensors  on 
DMSP  provide  a  rich  and  complex  source  for  observing 
electrodynamics  of  M-I  coupling  under  highly  driven 
conditions.  For  clarity  we  divide  the  observations  into 
four  parts.  The  first  subsection  contains  information  about 
interplanetary  drivers  and  the  global  geospace  response. 
The  latter  is  manifest  through  the  evolution  of  standard 
geomagnetic  indices  and  the  movement  of  auroral  bound¬ 
aries  [Newell  et  al.,  1996]  during  the  main  phase.  The 
second  subsection  presents  a  synopsis  of  detailed  particle 


Figure  2.  The  64-s  averaged  values  of  (top)  solar  wind 
density  (dotted  line)  and  speed  (solid  line),  (middle)  IMF 
By,  and  (bottom)  IMF  Bz  measured  by  the  ACE  satellite 
near  LI  on  6  April  2000. 

and  field  measurements  associated  with  intense  FAC  sheets. 
We  show  an  example  of  FAC  pairs  dividing  into  three  sheets 
as  they  relax  toward  quiet  levels.  The  third  subsection 
considers  circuit  impedances  estimated  from  observed  fluxes 
of  precipitating  electrons  [Robinson  et  al.,  1987]  and  ratios 
of  highly  correlated  Ey  and  6BZ  variations  in  equation  (7). 
The  Ey  and  6BZ  data  are  considered  from  the  viewpoint  of 
simple  DC  circuits  and  transmission  lines.  The  fourth 
subsection  concerns  ground  magnetometer  measurements 
taken  at  the  times  of  the  intense  FACs. 

4.1.  Geophysical  Background 

[29]  Partial  overviews  of  interplanetary  drivers  and  of 
geospace  responses  during  the  last  10  hours  of  6  April  2000 
are  provided  in  Figures  2  and  3,  respectively.  Figure  2 
shows  64  s  averaged  values  of  (top)  solar  wind  density  and 
speed,  (middle)  IMF  By,  and  (bottom)  IMF  Bz  as  measured 
by  plasma  and  magnetic  field  sensors  on  the  ACE  satellite 
near  the  first  Lagrangian  point  (Li).  These  data  show  that  at 
1604  UT  the  solar  wind  speed  increased  from  ~400  km/s  to 
steady  values  near  575  km/s.  At  the  same  time  the  density 
rose  from  6  to  20  cm-3.  Near  the  end  of  the  day,  the  density 
rose  to  ~40  cm-3.  After  1604  UT  the  IMF  had  strongly 
negative  Y  and  Z  components.  IMF  Bx  (not  shown)  was 
near  zero  all  day.  The  interplanetary  shock  passed  the  Wind 
satellite,  ~55  RE  upstream  of  the  Earth,  shortly  before 
1630  UT.  Compressional  waves  caused  by  the  shock’s 
impact  on  the  magnetopause  were  detected  by  ground 
magnetometers  at  ~1640  UT. 

[30]  Figure  3  shows  (top)  SYM-H  and  (bottom)  AE 
indices  for  6  April  2000  (available  at  http://swdcdb.kugi. 
kyoto-u.ac.jp/index.html).  SYM-H  summarizes  perturba¬ 
tions  in  the  horizontal  component  of  the  Earth’s  magnetic 
field  as  measured  around  the  world  at  midlatitude  to  low- 
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Figure  3.  One-minute  averaged  values  of  the  (top)  SYM-H 
and  (bottom)  auroral  electrojet  (AE)  indices. 


latitude  stations.  Information  provided  by  the  SYM-H  index 
is  similar  to  Dst  but  is  presented  at  a  1-min  rather  than  a 
1-hour  cadence.  Values  of  the  AE  index  are  also  presented 
as  1-min  averages.  SYM-H  rose  from  —27  to  +19  nT 
between  1 640  and  1 645  UT,  then  steadily  decreased  through 
the  rest  of  the  day  to  300  nT  as  the  ring  current 
intensified.  Sustained  recovery  began  early  on  7  April. 
The  auroral  electrojet  intensified  immediately  after  the 
sudden  storm  commencement  (SSC)  event.  Activity  was 
sporadic,  three  times  exceeding  2000  nT  before  declining 
after  1750  UT.  A  10-min  auroral  activation  began  at 
2000  UT,  with  AE  reaching  1085  nT.  This  excursion  was 
dominated  by  much  larger  contributions  from  AL  than  AU, 
suggesting  a  brief  episode  of  substorm-like  activity.  Rela¬ 
tive  auroral  quieting  with  AE  ranging  between  600  and 
200  nT  marked  the  remainder  of  the  day.  The  interval  of 
concern  extends  from  1900  to  2400  UT  when  SYM-H  was 
at  or  near  its  most  negative  values. 

[31]  Figure  4  shows  the  responses  of  auroral  boundaries 
[Newell  et  al.,  1996]  to  interplanetary  fluctuations.  Plotted 
are  the  magnetic  latitudes  of  Bl,  B2e,  B2i,  and  B5  observed 
by  five  DMSP  satellites  in  the  evening  and  dawn  MLT 
sectors  between  1900  and  2130  UT.  There  are  two  reasons 
for  the  lack  of  auroral  boundaries  on  the  morning  side  prior 
to  ~2020  UT.  First,  a  data  gap  in  DMSP  13  measurements 
occurred  between  1730  and  2000  UT.  Second,  the  other 
DMSP  satellites  crossed  the  auroral  oval  well  onto  the 
dayside  and  disallow  meaningful  applications  of  nightside 
boundaries  [Newell  et  al.,  1996],  Dashed  lines  on  both  plots 
represent  least  square  fits  to  the  magnetic  latitudes  of  the 
boundaries  as  functions  of  UT.  They  serve  only  as  guides 
for  the  eye.  The  data  illustrate  that  in  both  the  evening  and 
dawn  sectors,  the  magnetic  latitudes  of  Bl  and  B2e  were 
fairly  constant  during  the  period  of  interest.  B5  and  B2i 
systematically  shifted  toward  lower  latitudes,  indicating  that 
the  polar  cap  expanded  and  the  inner  edge  of  quasi-dipolar 
field  lines  moved  earthward  as  the  main  phase  progressed. 
Near  2020  UT,  B2i  locations  in  the  evening  sector  are  well 
equatorward  of  the  trend  line  at  MLat  <60°,  suggesting  that 
the  inner  magneto  tail  reconfigured  near  this  time. 

4.2.  Intense  Field-Aligned  Currents 

[32]  Figure  5  shows  particle  and  field  data  acquired  by 
(a)  F15  and  (b)  F12  in  the  evening  sector  and  (c)  FI 3  in  the 


dawn  sector,  respectively.  From  top  to  bottom  the 
panels  for  each  satellite  show  (1)  directional  energy  flux 
(E-Flux)  in  eV  cm-2  s-1  sr-1  of  down-coming  electrons 
(upper  line)  and  ions  (lower  line);  (2)  directional  number 
fluxes  (N-Flux)  of  down-coming  electrons  (upper  line) 
and  ions  (bottom  line)  with  energies  between  30  eV  and 
30  keV,  in  cm-2  s-1  sr-'1;  (3)  average  energy  of 
precipitating  ions  (upper  line)  and  electrons  (lower  line), 
defined  as  EAVe  =  E-Flux/N-Flux  (calculated  this  way, 
EAVe  is  about  twice  the  mean  thermal  energy  of  the 
species  population  [Robinson  et  al.,  1987]);  (4)  the  cross¬ 
track  component  of  the  magnetic  field  6BZ  in  nanotesla; 
and  (5)  the  in-track  component  of  the  electric  field  EY 
calculated  from  the  observed  ion  drift  and  IGRF  magnetic 
field  values  [Rich  and  Hairston,  1994],  Vertical  dashed 
lines  indicate  times  when  the  DMSP  satellites  crossed  Bl, 
B2e,  B2i,  and  B5  [Newell  et  al.,  1996].  Data  are  plotted 
as  functions  of  UT,  MLat,  and  MLT. 

[33]  Figure  5a  shows  that  FI 5  entered  the  evening  sector 
auroral  oval  from  the  southern  polar  cap  (B5)  near 
2015:30  UT.  A  rapid  rise  in  electron  fluxes  (panel  2)  and 
a  reversal  in  EY  (panel  5)  mark  the  transition.  Attention  is 
directed  to  the  ~1500  nT  decrease  in  6BZ  detected  between 
2018:10  and  2019:30  UT  (panel  4).  This  nearly  uniform 
decrease  in  6BZ  with  time  corresponds  to  a  2  pA/m2 
outward  FAC  or  a  net  electron  flux  of  ~109  cm-2  s-1  into 
the  ionosphere.  Figure  5a  also  shows  that  down-coming 
electron  fluxes  were  >109  cm-2  s-1  sr-1  in  the  upward  FAC 
sheet  (panel  2).  Assuming  isotropy  over  the  down-coming 
hemisphere,  we  see  that  electrons  with  energies  between 
30  eV  and  30  keV  can  account  for  the  observed  upward 
FAC.  The  EAVe  plot  (panel  3)  indicates  that  most  current- 
carrying  electrons  had  E  <  500  eV.  The  largest  6BZ 
deflection,  marking  the  interface  between  upward  and 
downward  FAC  sheets,  occurred  near  —56°  MLat  between 
the  B2i  and  B2e  boundaries.  A  deflection  of  1500  nT 
corresponds  to  an  integrated  FAC  J||  of  ~1.2  A/m  or 
~1.2  MA  for  every  1000  km  extent  in  longitude  (~1  hour 
in  LT).  Current  continuity  requires  that  a  height-integrated 
Pedersen  current  Jp,  equal  to  Jp,  flows  across  the  interface 
between  the  upward  and  downward  FAC  sheets. 

[34]  Near  the  transition  between  the  upward  and  down¬ 
ward  FAC  sheets,  a  maximum  electric  field  of  ~40  mV/m 
was  inferred  from  drift  meter  measurements  and  the  local 
IGRF  magnetic  field.  To  facilitate  direct  comparisons  be- 
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Figure  4.  Locations  of  Bl,  B2e,  B2i,  and  B5  auroral 
boundaries  observed  by  DMSP  satellites  in  the  (left) 
evening  and  (right)  morning  MLT  sectors  between  1900 
and  2200  UT  on  6  April  2000. 
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Figure  5.  Particle  and  field  measurements  by  (a)  DMSP  FI 5  and  (b)  FI 2  in  the  evening  sector  after 
2015  UT,  and  (c)  F13  in  the  morning  sector  after  2024  UT  on  6  April  2000.  From  top  to  bottom  the  panels 
for  each  satellite  show:  (1)  directional  energy  flux  (E-Flux)  in  eV  cm'2  s”1  sr'1  of  down-coming 
electrons  (upper  line)  and  ions  (lower  line);  (2)  directional  number  fluxes  (N-Flux)  of  down-coming 
electrons  (upper  line)  and  ions  (bottom  line)  with  energies  between  30  eV  and  30  keV,  in  cm-2  s_1  sr“'; 
(3)  average  energy  of  precipitating  ions  (upper  line)  and  electrons  (lower  line),  defined  as  E^ve  = 
E-Flux/N-Flux;  (4)  the  cross-track  component  of  the  magnetic  field  6BZ  in  nanotesla;  and  (5)  the  in-track 
component  of  the  electric  field  EY  calculated  from  the  observed  ion  drift  and  IGRF  magnetic  field  values 
[Rich  and  Hairston,  1994].  Vertical  dashed  lines  indicate  times  when  the  DMSP  satellites  crossed  Bl ,  B2e, 
B2i,  and  B5  [Newell  et  al.,  1996].  Data  are  plotted  as  functions  of  UT,  MLat,  and  MLT.  See  color  version 
of  this  figure  at  back  of  this  issue. 


tween  electric  fields  in  the  evening  and  dawn  MLT  sectors, 
we  have  plotted  EY  data  in  panel  5  of  Figures  5  a,  5b,  and  5c 
on  the  same  scale.  This  suppresses  information,  available  in 
standard  ion  drift  meter  displays,  that  indicate  a  1  km/s 
subauroral  ionospheric  drift  event  occurred  just  equatorward 
of  Bl  [Anderson  et  al.,  1993].  Magnetic  deflections  due  to 
R2  currents  and  subauroral  electric  fields  extend  equator- 
ward  to  approximately  -45°  MLat. 

[35]  Less  than  2  min  after  DMSP  FI  5,  F12  crossed  the 
auroral  oval  near  the  20  MLT  meridian,  where  it  detected 
similar  6BZ  and  EY  variations.  Particle  fluxes  measured  by 
the  SSJ4  sensor  on  F12,  shown  in  Figure  5b  panel  2,  yield 
one  additional  relevant  observation  that  differs  significantly 
from  the  FI 5  observations.  Data  plotted  in  Figure  5a 
panels  1  and  2  show  that  on  FI 5  in  the  downward  FAC 
region,  energy  fluxes  of  down-coming  ring  current  ions 
exceeded  those  of  precipitating  electrons.  Electron  number 
fluxes  were  always  greater  than  ion  fluxes.  DMSP  F12 
measurements  in  Figure  5b  panels  1  and  2  indicate  that 
across  the  downward  current  sheet  (MLat  >—56°)  ion 
number  fluxes  were  larger  than  those  of  electrons  by  a 
factor  of  two  or  more.  The  difference  is  attributed  to  intense 


fluxes  of  low  energy  (<1  keV)  ions  that  the  SSJ4  sensor  on 
FI 2  could  detect,  but  the  one  on  FI 5  could  not. 

[36]  Figure  6  shows  DMSP  F12  measurements  acquired 
between  2020  and  2023  UT  in  greater  detail.  Figure  6a 
repeats  the  6BZ  trace.  Figure  6b  shows  the  difference 
between  electron  and  ion  number  fluxes.  The  vertical  dashed 
line  marks  the  approximate  location  of  the  transition  from 
upward  to  downward  FACs.  To  the  left  of  this  line,  where  the 
6BZ  trace  has  a  negative  slope  indicating  a  current  out  of  the 
ionosphere,  the  flux  of  precipitating  electrons  exceeds  that  of 
the  ions  by  ~3  x  108/(cm2  s  sr).  Near  2022  UT  the  opposite 
FAC  polarity  is  observed  and  the  precipitating  ion  flux 
dominates.  Figures  6c  and  6d  show  5-s  averaged  electron 
and  ion  spectra  centered  near  2021:00  and  2021:51  UT  in 
the  upward  and  downward  FAC  sheets,  respectively.  In  the 
upward  FAC  sheet,  low  (<1  keV)  electrons  dominate  the 
spectrum.  In  the  downward  FAC  sheet  the  flux  of  <1  keV 
ions  almost  matches  that  of  the  electrons.  At  energies 
>2  keV  the  ion  spectrum  is  relatively  flat,  indicating  the 
presence  of  substantial  fluxes  at  energies  >30  keV,  the 
upper  limit  of  SSJ4  detection.  As  mentioned  above,  DMSP 
satellites  cannot  monitor  upward  electron  or  ion  fluxes. 
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Figure  6.  Particle  and  field  characteristics  of  the  FAC 
sheets  crossed  by  DMSP  F12  between  2020  and  2023  UT 
on  6  April  2000.  Panels  show  (a)  variations  of  8BZ  and 
(b)  differences  between  electron  and  ion  number  fluxes  as  a 
function  of  UT,  as  well  as  5-s  averaged  ion  and  electron 
spectra  measured  by  SSJ4  within  the  (c)  upward  and 
(d)  downward  FAC  sheets.  See  color  version  of  this  figure 
at  back  of  this  issue. 

However,  the  DMSP  F12  observations  in  Figures  6b  and  6d 
strongly  suggest  that  precipitating  ions  dominated  over 
precipitating  electrons. 

[37]  Returning  to  Figure  5c,  we  see  that  DMSP  FI 3  data 
were  acquired  near  the  dawn  meridian  in  the  northern 
ionosphere.  Electron  fluxes  (panel  2)  rose  above  polar 
rain  levels  at  ~2025:50  UT.  The  polarity  of  Ey  reversed 
~30  s  later.  F13  entered  the  main  oval  (B5)  just  before 
2027:00  UT.  In  the  region  poleward  of  B5,  electron  and  ion 
fluxes  were  often  near  background  levels,  leading  to  large, 
nongeophysical  variations  in  calculated  values  of  EAVe- 
Since  the  polar  cap  is  outside  the  domain  of  concern  for 
this  paper,  we  do  not  display  EAVe  values  in  this  region.  A 
comparison  with  data  in  Figure  5a  panels  1  and  2  indicates 
that  BPS  electron  precipitation  was  more  structured  in  the 
dawn  than  in  the  evening  sector  [Winningham  et  al.,  1975]. 
Particle  fluxes  observed  in  the  morning  sector  by  Fll  and 
F14  (not  shown)  were  similar  to  those  of  F13.  Electron 
energy  fluxes  >1  erg  cm-2  s_1  sr~‘  were  observed  for 
more  than  a  minute  after  2027  and  2030  UT.  A  pair  of 
intense  FAC  sheets,  with  a  magnetic  deflection  >1400  nT 
(~1 . 1  A/m),  was  observed  near  2030  UT  (panel  4).  A  sharp 
intensification  of  electron  fluxes  (panels  1  and  2)  marked 
the  transition  from  downward  to  upward  FAC.  In  this  case, 
Ey  «  90  mV/m  was  inferred  from  drift  meter  measurements 
near  the  maximum  6BZ.  In  the  region  of  the  downward 
FAC,  low  ion  fluxes  indicate  that  upgoing  ionospheric 
electrons  were  probably  the  prime  current  carriers.  Again 
the  maximum  deflections  in  Ey  and  6BZ  occurred  between 
B2e  and  B2i.  Finally,  we  note  that  in  both  the  morning  and 
evening  sectors,  R2  FACs  extended  to  subauroral  latitudes. 

[38]  An  intriguing  aspect  of  the  satellite  observations 
concerns  8BZ  signatures  as  the  FAC  system  relaxed. 
Figure  7  shows  8BZ  signatures  observed  by  (top)  F15  and 
(bottom)  FI 2  during  the  next  evening  pass  over  the  northern 
auroral  oval.  As  the  two  spacecraft  moved  poleward  into 


and  across  the  evening  sector  of  the  northern  auroral  oval, 
they  detected  three,  rather  than  two,  current  sheets,  as 
indicated  by  the  arrows  in  Figure  7.  At  low  magnetic 
latitudes  both  satellites  encountered  FACs  directed  out  of 
the  ionosphere.  Poleward  of  this  upward  current  sheet  they 
crossed  the  standard,  evening-sector  R1/R2  system.  No  sign 
of  an  additional  FAC  sheet  appeared  in  FI 3  and  F14 
measurements  as  the  satellites  crossed  the  southern  auroral 
oval  in  the  morning  sector.  There  are  five  empirical  facts  to 
be  noted.  (1)  The  equatorward  FACs  out  of  the  ionosphere 
occurred  between  B1  and  B2e.  (2)  DMSP  FI 5  and  F12 
measured  average  upward  current  densities  of  ~  1.4  pA/m2 
and  1 .0  pA/m2,  with  integrated  currents  of  0.56  Am  and 
0.4  Am,  respectively.  (3)  Within  the  equatorward  upward 
current  sheet,  electric  fields  were  weak  and  equatorward. 
(4)  Maximum  6BZ  deflections  at  the  R1/R2  interface  oc¬ 
curred  between  B2e  and  B2i.  (5)  Total  6BZ  perturbations 
measured  by  FI 5  and  F12  in  the  R2  currents  were  ~1400  nT 
(1.1  Am)  and  1000  nT  (0.8  Am),  respectively. 

[39]  Particles  monitored  by  FI 2  and  FI 5  between  B1 
and  B2e  for  this  pass  (not  displayed)  indicate  that  electron 
N-fluxes  always  exceeded  those  of  ions.  In  the  previous 
southern  hemisphere  pass,  FI 2  measurements  show  that  the 
opposite  relationship  prevailed  (Figure  6).  However,  ion  E- 
fluxes  from  the  central  plasma  sheet  (between  B2e  and  B2i) 
were  similar  to  those  observed  during  the  southern  auroral 
oval  crossing.  The  discussion  section  considers  implications 
of  these  measurements  for  our  understanding  of  ring  current 
pressure  distributions. 

[40]  Analysis  of  DMSP  data  acquired  during  the  magnetic 
storm  was  facilitated  by  the  relative  spacecraft  orbital 
phasing  prevailing  on  6  April  2000.  Fll,  FI 3,  and  F14 
crossed  the  high-latitude  ionosphere  in  the  same  hemisphere 
within  2  min  of  each  other.  The  F12  and  FI  5  spacecraft  flew 
over  the  conjugate  polar  ionosphere  ~10  min  earlier  than 
Fll,  F13,  and  F14  but  within  1.5  min  of  each  other.  Table  1 
lists  the  UT  as  well  as  the  magnetic  and  geographic 
coordinates  at  which  the  DMSP  satellites  crossed  the  largest 
FAC  perturbations.  Value  of  the  largest  amplitude  of  |8BZ| 


Ml  H2e  I12i  U5 


II  r  2(1:48  20:51  20:54 

MUt  47.0  58.5  69.6 

MU'  21.1  21.1  21.2 


HI  112c  I12i 


Ml.al  35.2  49.7  61.2 

Mil  20.5  20.5  20.5 

Figure  7.  Values  of  SBZ  observed  by  (top)  FI 5  and 
(bottom)  F12  while  crossing  the  northern  auroral  oval 
between  2048  and  2054  UT.  Upward/downward  pointing 
arrows  indicate  the  polarities  of  FACs  out  of/into  the 
ionosphere. 
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Table  1.  Extrema  of  6BZ  and  EY 


DMSP 

UT 

MLAT 

MLT 

MLONG 

GLAT 

GLONG 

|5BZ|,  nT 

|EvL  mV/m 

F12 

2012.2 

-56.0° 

19.9 

66° 

-55.0° 

14.5° 

1326 

41.5 

F13 

2030.3 

58.7° 

6.7 

225° 

-63.0° 

159.9° 

1424 

91.5 

F14 

2030.7 

59.0° 

8.4 

255° 

61.2° 

199.8° 

1354 

82.9 

F15 

2019.4 

-55.9° 

20.9 

79° 

-52.1° 

26.4° 

1506 

42.5 

appears  in  the  eighth  column.  Maxima  for  |EY|  were 
acquired  spatially  close  to,  but  not  necessarily  coincident 
with,  6BZ  extrema.  Equation  (6)  indicates  that  separations  in 
|8BZ|  and  |EY|  extrema  reflect  the  presence  of  conductivity 
gradients  near  the  upward/downward  FAC  interface. 

[41]  Three  empirical  aspects  of  data  in  Table  1  call  for 
comment.  First,  8BZ  extrema  measured  by  all  of  the  satellites 
corresponded  to  FACs  with  Jm  >  1  A/m.  Second,  all  8BZ 
extrema  were  observed  at  magnetic  latitudes  equatorward  of 
60°.  Third,  electric  fields  detected  at  similar  MLTs  had 
comparable  values.  The  latitudinal  distributions  Of  Ey 
detected  by  Fll  and  F13  were  nearly  identical  across  the 
morning  auroral  oval.  Note,  however,  that  EY  values  mea¬ 
sured  in  the  evening  sector  were  about  half  those  detected 
on  the  morning  side  (compare  panel  5  of  Figures  5a,  5b,  and 
5c).  Ohm’s  law  dictates  that  ionospheric  conductance  was 
larger  in  the  evening  than  the  morning  sector  by  about  a 
factor  of  two. 

[42]  In  the  remaining  analysis  of  DMSP  data  we  adopt 
|6Bz(max  as  our  organizing  parameter  for  understanding  M-I 
coupling.  This  choice  is  dictated  by  two  considerations. 
First,  equation  (5)  shows  that  |6Bz|max  is  directly  propor¬ 
tional  to  J||  the  integrated  FAC  coupling  the  magnetosphere 
to  the  ionosphere.  As  such  it  maps  to  the  ionospheric 
footprint  of  maximum  Pedersen  current.  Second,  during  the 
magnetic  storm’s  main  phase,  |8Bz|max  was  observed  at 
relatively  low  magnetic  latitudes  that  normally  map  to  the 
inner  magnetosphere,  which  primarily  acts  as  a  current 
generator  [Vasyliunas,  1970].  Figure  8  plots  all  57  detections 
of  |6Bz|max  observed  by  the  four  DMSP  satellites  carrying 
SSM  between  1600  and  2400  UT.  For  reference,  Figure  8 
also  contains  a  plot  of  the  AE  index.  Somewhat  arbitrarily 
we  refer  to  FAC  events  as  intense  if  |8Bz|max  >  1000  nT.  We 
see  that  of  the  14  intense  events,  11  occurred  late  in  the 
main  phase  after  2215  UT  when  AE  was  generally 
decreasing  after  several  brief  intensifications  and  SYM-H 
was  either  decreasing  or  oscillating  near  its  minimum  value 
(Figure  3). 

[43]  Vertical  lines  on  Figure  8  at  2015,  2155,  2250,  and 
2345  UT  highlight  intensifications  of  the  FAC  system.  The 
events  described  in  Figure  5  and  in  Table  1  occurred  near 
2015  UT,  shortly  after  AE  decreased  by  ~500  nT 
(Figure  2).  The  FAC  intensifications  appear  to  occur  on 
fairly  rapid  time  scales.  As  F13  and  F14  crossed  the 
northern  auroral  oval  they  detected  magnetic  perturbations 
of  646  nT  at  2012  UT  (MLat  =  57.7°,  MLT  =  18.2)  and 
797  nT  at  2014  UT  (MLat  =  59.9°,  MLT  =  20.6),  respec¬ 
tively.  At  2019  UT,  F15  observed  a  1506  nT  perturbation  in 
the  conjugate  ionosphere  (MLat  =  -55.9°,  MLT  =  20.9); 
2  min  later  FI 2  also  crossed  the  same  intense  FAC  sheet.  A 
similar  rapid  intensification  of  FACs  occurred  shortly 
after  2200  UT  when  three  satellites  measured  |8Bz|max  > 
1200  nT.  Note,  however,  that  at  ~220 1  UT,  F12  and  F15 
detected  maximum  perturbations  of  769  and  1317  nT, 
respectively,  while  separated  by  <2°  in  MLat  and  ~1  hour 


in  MLT.  While  intensifications  of  FACs  appear  four  times 
during  the  interval  plotted  in  Figure  8,  no  evidence  of 
sawtooth  oscillations  appears  in  particle  fluxes  observed  at 
geosynchronous  orbit  [Huang  et  al.,  2003]. 

4.3.  Circuit  Conductance 

[44]  The  previous  section  derives  relationships  between 
variations  of  6BZ  and  Ey  required  to  maintain  current 
continuity  in  the  ionosphere  beneath  DMSP  satellites. 
Smiddy  et  al.  [1980]  showed  that  in  regions  of  nearly 
uniform  conductance  equation  (6)  reduces  to  EP  =  A8BZ/ 
[poA  Ey],  The  symbol  A  indicates  changes  in  8BZ  or  Ey.  In 
practical  units  EP  (mho)  =  A6BZ  (nT)/[1.256  •  AEy  (mV/ 
m)].  These  equations  assume  that  the  infinite  current  sheet 
approximation  is  valid.  For  the  events  of  Figure  5,  simul¬ 
taneous  6BZ  and  5By  variations  were  nearly  monotonic  and 
highly  correlated,  indicating  that  the  spacecraft  indeed 
crossed  quasi-infinite  current  sheets  at  nearly  constant 
attack  angles  ct  =  Tan-'  (6Bz/8By).  Values  of  a  listed  in 
the  last  column  of  Table  2  show  that  the  DMSP  satellites 
crossed  FAC  sheets  at  nearly  normal  incidence.  Fung  and 
Hoffman  [1992]  addressed  the  more  general  case  where  end 
effects  or  filamentary  structures  are  important. 

[45]  Table  2  also  summarizes  information  contained  in 
Figure  5.  Listed  data  indicate  the  beginning  UTj  and 
duration  At  of  intervals  in  which  6BZ  and  Ey  variations 
were  highly  correlated,  consistent  with  crossing  regions  of 
nearly  uniform  EP.  Correlation  coefficients  within  the  cur¬ 
rent  sheets  were  >0.98.  Values  of  EP  estimated  from  8BZ 
and  Ey  variations  within  the  intense  upward  and  downward 
FAC  sheets  are  listed  in  column  8  of  Table  2.  Except  in  the 
downward  FAC  sheet  crossed  by  F13,  estimated  values  of 
EP  were  >25  mho. 

[46]  Measurements  from  the  SSJ4  sensor  allow  indepen¬ 
dent  estimates  of  the  ionospheric  conductance  using  empir¬ 
ical  formulas  suggested  by  Robinson  et  al.  [1987]. 

EP  =  [40  Eave  -\/($E)]/[l6  +  Eave]  (12)' 

and 

Eh  =  0.45  Ep(EAve)0  85.  (13) 


Figure  8.  Values  of  AE  (solid  line)  and  |SBz|max  (circles) 
plotted  as  a  function  of  UT. 
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Table  2.  Auroral  FAC  Circuit  Parameters 


DMSP 

UT, 

At,  s 

6Bzr,  nT 

6BZF,  nT 

Eyi,  mV/m 

Eyf,  mV/m 

Ep,  mho 

a 

F13 

2028:46 

64 

600 

1340 

3.0 

105.4 

5.9 

68.3° 

F13 

2030:30 

21 

1370 

498 

48.5 

25.4 

30.0 

83.9° 

F15 

2018:22 

64 

-218 

-1487 

0.4 

37.4 

27.3 

84.5° 

F15 

2019:30 

76 

-1469 

-496 

34.7 

9.8 

31.1 

89.2° 

Here  EAVe  represents  the  average  energy  of  precipitating 
electrons  in  keV  and  $e  is  the  electron  energy  flux  in  ergs/ 
(cm2  s).  To  help  estimate  ionospheric  conductances  we 
assume  that  electron  distribution  functions  were  isotropic 
over  the  down-coming  hemisphere.  Calculated  values  of  £P 
(top)  and  EH  (bottom)  due  to  energetic  electron  precipita¬ 
tion  observed  by  F15  (left)  and  FI 3  (right)  are  plotted  as 
functions  of  UT  in  Figure  9.  The  vertical  dashed  lines  mark 
the  intervals  during  which  £P  was  estimated  from  highly 
correlated  6BZ  and  EY  variations.  The  downward  and 
upward  pointing  arrows  located  between  the  dashed  lines 
signify  the  directions  of  FACs  into  and  out  of  the 
ionosphere,  respectively.  Heavy  horizontal  bars  within  the 
dashed  lines  indicate  calculated  values  of  £P  listed  in 
Table  2.  Comparing  Pedersen  conductances  estimated  using 
two  different  techniques,  we  see  that  the  average  values  of 
£P  estimated  from  equation  (12)  were  significantly  lower 
than  those  obtained  from  equation  (7).  Attention  is  directed 
to  the  upward  current  sheets  where  observed  variations  of 
6BZ  and  EY  were  monotonic  (Figure  5)  and  highly 
correlated.  Here  equation  (6)  indicates  that  £P  should  have 
been  nearly  constant.  In  both  instances  £P  plots  estimated 
from  electron  measurements  suggest  the  presence  of 
gradients.  We  defer  comment  on  these  discrepancies  to  the 
discussion  section. 

4.4.  Ground  Measurements 

[47]  Black  circles  on  the  trajectories  of  Figure  1  and 
spacecraft  designation  in  Figure  10  indicate  the  locations 
of  |6BZ|  and  |EY|  extrema  in  magnetic  and  geographical 
coordinates.  Table  1  lists  the  exact  locations.  F12  and  FI 5 
were  flying  over  the  Indian  Ocean  to  the  southeast  of  the 
African  continent.  F13  and  F14  were  over  northeastern 
Siberia.  The  8BZ  extrema  map  to  the  ionospheric  location 
of  the  most  intense  Pedersen  currents  linking  downward  and 
upward  FAC  sheets.  A  priori,  one  might  expect  that  mag¬ 
netometers  at  ground  stations  located  beneath  these  currents 
would  detect  comparably  large  magnetic  fluctuations.  Cape 
Schmidt  (CPS)  at  68.9°N,  180.5°E  provides  the  only 
magnetic  field  measurements  presently  available  from  the 
Siberian  sector.  The  station  lies  about  half  way  between  the 
longitudes  of  F13  and  FI 4.  However,  it  is  ~6°  in  latitude  to 
the  north  of  these  FAC  encounters.  CPS  records  show  only 
modest  (~200  nT)  variations  between  2010  and  2100  UT. 

[48]  The  locations  of  F12  and  F15  over  the  Indian  Ocean 
make  direct  comparisons  with  subsatellite  magnetometers 
impossible.  No  significant  magnetic  perturbations  were 
recorded  at  stations  arrayed  along  the  coast  of  Antarctica. 
Table  1  gives  the  magnetic  latitudes  and  longitudes  of  F12 
(-56°,  66°)  and  F15  (-55.9°,  19°).  If  the  inner  magneto¬ 
sphere  acts  as  a  current  generator  [Kan  and  Sun,  1985], 
spacecraft  flying  above  the  conjugate  ionosphere  should 
observe  8BZ  signatures  similar  to  those  detected  by  F12  and 
FI 5.  Consequently,  similar  magnetic  perturbations  would 


be  detected  at  conjugate  ground  locations.  The  magnetic 
latitude  and  longitude  (56.3°,  74.7°)  of  the  Valentia  (VAL) 
station  in  Ireland  is  nearly  conjugate  to  the  places  where 
F12  and  FI 3  detected  6BZ  extrema.  Figure  11  shows  the 
magnitude  and  three  components  of  magnetic  perturbations 
observed  at  VAL  between  1400  and  2100  UT.  Variations 
measured  between  2000  and  2100  UT  were  <100  nT.  An 
examination  of  records  from  a  large  number  of  auroral  and 
middle  latitude  stations  shows  that  magnetic  variations 
commensurate  with  those  observed  by  DMSP  were  not 
detected  on  the  ground.  Interested  readers  may  verify  this 
conclusion  by  examining  the  many  ground  magnetic 
records  for  6  April  2000,  available  on  the  Kyoto  web  site 
(http://swdcdb.kugi.kyoto-u.ac.jp/index.html). 

5.  Discussion 

[49]  The  previous  section  shows  particle  and  field  mea¬ 
surements  acquired  during  the  main  phase  of  the  6  April 
2000  magnetic  storm,  focusing  on  nearly  simultaneous 
observations  of  intense  (Jy  >  1  A/m)  FACs  in  the  evening 
and  dawn  MLT  sectors.  Large  magnetic  perturbations, 
marking  the  R1/R2  transition,  occurred  between  the  B2e 
and  B2i  auroral  boundaries  [Newell  et  al.,  1996].  Height- 
integrated  conductivities  estimated  from  electric  and 
magnetic  field  variations  [Smiddy  et  al.,  1980]  were 
systematically  larger  than  those  determined  from  electron 
fluxes  [Robinson  et  al.,  1987].  Larger  electric  fields 
appeared  in  the  morning  than  in  the  evening  sector,  with 
similar  FAC  inputs.  Finally,  available  data  from  ground 
stations  show  no  commensurate  magnetic  perturbations 
while  DMSP  was  detecting  intense  FACs. 


Figure  9.  (Top)  values  of  EP  and  (bottom)  £H  estimated 
from  energetic  electron  fluxes  measured  by  (left)  FI 5  and 
(right)  FI  3  plotted  as  functions  of  UT.  The  heavy  horizontal 
lines  indicate  values  of  EP  estimated  from  6BZ  and  EY 
variations. 
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Distribution  of  Geomagnetic  Observatories 


Figure  10.  World  map  showing  the  positions  of  four  DMSP  satellites  when  they  encountered  extrema 
of  |8BZ|  and  |Ey|.  Locations  of  the  Cape  Schmidt  (CPS),  Memambetsu  (MEM),  and  Valentia  (VAL) 
magnetic  stations  are  indicated  for  reference. 


[50]  Observations  presented  here  are  unusual  but  not 
unique.  An  informal  survey  of  DMSP  measurements 
revealed  that  J||  >  1  A/m  events  regularly  occur  during  the 
main  phases  of  magnetic  storms  with  Dst  <  —200  nT. 
Distinct  ground  signatures  accompanied  many,  but  not  all, 
of  the  intense  stormtime  FAC  events.  At  the  suggestion  of 
G.  Lu  (private  communication,  2002),  we  examined  DMSP 
particle  and  field  data  acquired  during  the  superstorms  of 
15  July  2000  and  31  March  2001.  DMSP  measurements 
during  both  storms  bore  striking  similarities  to  the 
observations  reported  here.  Again  intense  FACs  observed 
at  840  km  were  unmatched  by  strong  magnetic  perturba¬ 
tions  at  the  ground.  In  both  of  these  superstorms  the  intense 
FACs  were  observed  at  latitudes  well  equatorward  of  the 
normal  auroral  oval,  and  EAVe  <  500  eV  for  precipitating 
electrons. 

[51]  This  section  discusses  (1)  some  implications  of 
DMSP  observations  for  understanding  the  dynamics  of 
magnetospheric  plasma  pressure  distributions,  (2)  applica¬ 
tions  of  Ohm’s  law  in  the  ionosphere,  (3)  the  absence  of  FAC 
magnetic  signatures  on  the  ground,  and  (4)  current  carriers. 

5.1.  Magnetospheric  Pressure  Distributions 

[52]  The  intense  FACs  crossed  by  DMSP  satellites  in  the 
evening  and  dawn  MLT  sectors  developed  during  the 
storm’s  main  phase.  Activity  detected  at  the  latitudes  of 
auroral  electrojet  stations  was  generally  quieting  but  was 
also  subject  to  sporadic  increase/decrease  cycles  (Figures  3 
and  8).  The  negative  slope  of  the  SYM-H  index  during  the 
period  of  interest  indicates  that  the  ring  current  intensified 
and/or  moved  closer  to  Earth.  The  fact  that  the  latitude  of 
B1  was  stable  and  that  5BZ  and  EY  were  confined  to  MLat 
>45°  suggests  that  the  ring  current’s  inner  edge  did  not 
move  much  further  earthward.  Moving  the  ring  current 
earthward  of  L  =  2  requires  that  a  significant  dawn-dusk 
electric  field  penetrate  more  closely  to  Earth.  This  inference 
is  consistent  with  CRRES  observations  that  ring  current 
ions  did  not  move  to  L  <  2  during  the  magnetic  storms  of 
June  and  July  1991  [ Burke  etal.,  1998;  Wilson  etal.,  2001], 


[53]  Vasyliunas  [1970]  showed  that  field-aligned  currents 
jll  flow  into  or  out  of  the  ionosphere  whenever  gradients  of 
the  magnetospheric  plasma  pressure  p  are  not  aligned  with 
gradients  in  magnetic  flux  tube  volumes  £,  =  J (ds/B) 


in 


—  (Bm/(2Be)be 


Vp  x  V 


(14) 


where  be  is  a  unit  vector  along  the  Earth’s  magnetic  field  at 
the  equator  (Be)  and  Bm  is  the  magnetic  field  at  the  top  of 
the  ionosphere.  V£,  has  a  strong  radially  outward  component 
at  all  local  times.  Magnetospheric  pressure  distributions 
estimated  by  Wing  and  Newell  [1998]  from  ion  fluxes 
observed  between  B2i  and  B5  indicate  that  pressure  maxima 
occur  in  the  midnight  to  late  evening  MLT  sector  at 
distances  from  the  Earth  that  depend  on  the  level  of 
geomagnetic  activity.  In  the  inner  magnetosphere  near  the 


Figure  11.  Three  components  and  intensity  of  magnetic 
perturbations  measured  at  Valentia  between  1400  and 
2100  UT  on  6  April  2000.  Numbers  listed  to  the  right  of 
the  plots  indicate  baseline  values  from  which  perturbations 
are  measured. 
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dusk  (dawn)  meridians  where  V£,  is  mostly  radial,  we 
expect  FACs  into  (out  of)  the  ionosphere  in  the  evening 
(morning)  sector  due  to  the  eastward  (westward)  pressure 
gradients. 

[54]  Beyond  the  B2i  boundary,  azimuthal  pressure  gra¬ 
dients  appear  weak  [Wing  and  Newell,  1998]  but  still 
maintain  maxima  close  to  the  midnight  meridian.  The 
location  of  the  largest  6BZ  and  EY  perturbations  between 
B2e  and  B2i  indicates  that  the  transition  between  R2  and  R1 
currents  occurred  in  the  inner  part  of  the  central  plasma 
sheet.  This  conclusion  agrees  with  the  work  of  Wing  and 
Newell  [2000]  who  show  that  R1  FACs  are  generated 
tailward  of  B2i  under  quiet  conditions  and  that  the  transition 
between  R1  and  R2  occurs  between  B2i  and  B2e. 

[55]  The  concept  of  magnetospheric  plasma  pressure 
gradients  controlling  the  distribution  of  ionospheric  FAC 
distributions  [Vasyliunas,  1970]  provides  guidance  for  un¬ 
derstanding  the  FAC  sheet  out  of  the  ionosphere  observed 
by  F12  and  F15  at  low  (~50°)  magnetic  latitudes  near 
2050  UT.  Comparisons  of  5BZ  measurements  in  the  south¬ 
ern  and  northern  hemispheres  (Figures  5  and  7)  show  that 
while  B2i  and  B5  moved  equatorward  by  ~2°,  the  magnetic 
latitudes  of  B1  and  B2e  were  steady.  The  ~650  nT  decrease 
in  8BZ  seen  by  F15  between  2049:00  and  2049:50  UT 
indicates  that  a  FAC  of  ~1.6  p,A/m2  flowed  out  of  the 
ionosphere.  The  average  intensity  of  the  downward 
FAC  encountered  between  2049:50  and  2050:50  UT  was 
~3  pA/m2,  about  twice  its  southern  hemisphere  value.  The 
outward  FAC  observed  by  F12  and  FI 5  between  B1  and 
B2e  requires  that  close  to  the  inner  edge  of  the  plasma  sheet 
near  2100  MLT,  Vp  had  a  westward  component.  In  the 
central  plasma  sheet,  between  the  mappings  of  B2e  and  B2i, 
the  eastward  component  of  Vp  intensified. 

[56]  The  requirement  for  a  westward  pressure  gradient  in 
the  innermost  plasma  sheet  suggests  the  following  scenario. 
The  brief  burst  of  substorm-like  activity  observed  near 
2000  UT  (cf.  AE  trace  in  Figure  7)  injected  plasma  close 
to  the  Earth  in  the  midnight  sector.  Because  of  the  short  time 
allowed  for  the  injection  and  ring  current  shielding  of 
magnetospheric  electric  fields,  the  pressure  maximum  was 
between  the  mapping  of  B2e  and  B2i.  Flowever,  a  large 
number  of  energetic  ions  reached  the  mapped  location  of  B 1 . 
When  the  substorm  driver  suddenly  turned  off,  gradient- 
curvature  drifts  dominated  ion  motion.  Most  of  the  pressure 
was  carried  by  westward  drifting  energetic  ions.  The  gradient 
drift  velocity  VG  of  an  energetic  ion  is 

VG  =  (pi/qB2)B  x  VB,  (15) 

where  p*  the  ion’s  magnetic  moment  is  a  constant  of  the 
motion.  In  a  magnetic  dipole  approximation  B  is  propor¬ 
tional  to  R-3,  where  R  is  the  distance  from  the  center  of  the 
Earth.  Thus  VG  is  proportional  to  R_1.  Since  VG  =  R  dip/dt, 
the  rate  at  which  energetic  ions  drift  in  local  time  to  the  west 
is  proportional  to  R”2.  Although  we  lack  information  about 
the  magnetic  field  structure  in  the  equatorial  plane,  our 
approximation  suggests  that  the  earthward  part  of  the 
injected  plasma  drifts  to  the  west  at  a  much  faster  angular 
speed  than  the  part  initially  further  from  the  Earth  in  the 
main  plasma  sheet  (poleward  of  B2e). 

[57]  If  dynamics  observed  during  the  magnetic  storm  of 
June  1991  were  repeated  on  6  April  2000,  then  inflationary 


effects  of  the  ring  current  should  have  caused  the  equatorial 
mapping  of  B1  to  lie  closer  to  6  than  to  3  RE  from  the  center 
of  the  Earth  [Burke  et  al.,  1998],  Under  these  circumstances, 
VB  should  have  been  unusually  strong  between  the  map¬ 
pings  of  B1  and  B2e.  Assuming  that  the  part  of  the  ring 
current  population  B 1  and  B2e  covered  a  relatively  narrow 
swath  in  local  time,  it  appears  plausible  that  the  centroid  of 
plasma  pressure  gradient-curvature  drifted  to  the  west  of 
the  DMSP  F12  and  FI 5  orbits.  In  the  evening  sector  Vp 
would  have  had  a  westward  component  to  drive  an  outward 
FAC.  From  the  perspective  of  current-carrying  particles, 
between  B1  and  B2e  SSJ4  detected  larger  (smaller)  fluxes 
of  precipitating  electrons  than  down-coming  ions  in  the 
northern  (southern)  hemisphere  pass. 

5.2.  Ionospheric  Ohm’s  Law 

[58]  Plots  in  Figure  9  compare  values  of  EP  calculated 
using  electron  fluxes  measured  by  the  SSJ4  sensors  on 
DMSP  FI 5  and  FI 3  [. Robinson  et  al.,  1987]  with  estimates 
based  on  simultaneous  electric  and  magnetic  field  variations 
observed  with  SSIES  and  SSM  [Smiddy  et  al.,  1980]. 
Except  in  the  region  of  upward  FACs  detected  by  F13, 
values  of  EP  estimated  using  the  electromagnetic  technique 
are  much  higher.  We  suggest  three  reasons  for  this  discrep¬ 
ancy.  First,  electromagnetic  calculations  of  EP  are  most 
valid  along  trajectory  segments  where  the  conductance  is 
nearly  constant.  Distributions  of  EP  determined  from  elec¬ 
tron  fluxes  suggest  that  this  condition  was  not  always  met. 
Second,  during  the  main  phase  of  this  storm  M-I  coupling 
conditions  exceeded  the  limits  of  validity  for  applying  the 
formulas  of  Robinson  et  al.  [1987].  Third,  the  equation 
derived  by  Smiddy  and  coworkers  assumes  that  the  M-I 
circuit  has  reached  a  steady  state.  F13  and  F14  detected 
relatively  small  values  of  |6Bz|max  in  the  evening  auroral 
oval  a  few  minutes  before  F12  and  F15  (Figure  8).  These 
observations  suggest  that  FI 5  data  presented  in  Figure  5a 
were  acquired  soon  after  the  intense  FAC  circuit  was 
activated.  It  takes  several  Alfven  bounce  periods  for  the 
M-I  circuit  to  reach  equilibrium  [Kan  and  Sun,  1985];  thus 
we  must  also  consider  transmission  line  effects.  We  address 
these  possibilities  in  turn. 

5.2.1.  Conductance  Gradients 

[59]  A  cursory  look  at  variations  in  energetic  electron 
fluxes  measured  by  DMSP  satellites  while  crossing  individ¬ 
ual  FAC  sheets  shows  that  EP  gradients  must  exist  on  small 
spatial  scale  sizes.  Indeed,  small-scale  variations  in  6BZ  and 
Ey  traces  attest  to  the  presence  of  similar-scale  EP  gradients. 
For  two  reasons  we  believe  that  conductance  gradients  are 
not  as  severe  as  suggested  by  estimates  presented  in  Figure  9. 
First,  in  choosing  regions  to  estimate  EP  from  6BZ  and  Ey 
variations,  we  required  that  the  correlation  between  the  two 
data  sets  exceed  0.98.  Second,  in  the  vicinity  of  6BZ  and  EY 
extrema,  conductances  calculated  using  the  formulas  of 
Robinson  et  al.  [1987]  are  inconsistent  with  simple  applica¬ 
tions  of  Ohm’s  law.  At  6BZ  extrema  J||  =  Ip  =  EPEY, 
independent  of  conductivity  profiles.  Here  EP  represents  the 
local  value  of  the  Pedersen  conductance  at  the  transition 
between  upward  and  downward  FACs.  In  these  regions  the 
formulas  of  Robinson  et  al.  [1987]  predict  EP  values  of  a 
few  mhos.  For  the  cases  shown  in  Figures  5a  and  5c,  Jji  « 
1.2  A/m.  F13  and  F15  measured  EY  values  of  91.5  and 
42.5  mV/m  (Table  1),  respectively.  Corresponding  EP  values 
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Table  3.  Alfven  Wave  Reflection  Coefficients 

DMSP/j ||  ASBz/AEy  SP/SA  R  EV,  mV/m  6B'Z|  nT 


F13  If  7.40  x  10-6  5.5  -0.70  338  435 

F13  H  3.79  x  10“5  28.4  -0.93  329  452 

F15  i[  3.43  x  1(T5  18.9  -0.90  369  668 


of  13.1  and  28.2  mho  exceed  those  shown  in  Figure  9 
between  the  upward  and  downward  current  sheets.  Thus,  EP 
gradients  alone  cannot  account  for  the  difference  between 
the  two  methods  of  estimation. 

5.2.2.  Applicability  Limits  of  Robinson  Formulas 

[60]  At  large  solar  zenith  angles  Robinson’s  formulas 
estimate  EP  and  EH  for  a  given  energy  flux  and  average 
energy  of  precipitating  electrons.  Robinson  et  al.  [1987] 
excluded  contributions  to  conductances  of  precipitating 
electrons  with  energies  <500  eV  and  all  protons  and 
cautions  against  use  of  this  formula  when  electron  energies 
are  low.  During  the  April  2000  magnetic  storm,  neglect  of 
these  conditions  leads  to  serious  underestimates  of  EP.  Data 
in  Figure  5  indicate  that  the  average  electron  energy  was 
below  500  eV.  These  electrons  are  stopped  through  colli¬ 
sions  with  thermospheric  neutrals  at  altitudes  above  190  km 
[Rees,  1964],  Ionization  created  by  these  incident  electrons 
places  the  Pedersen  current-closure  layer  at  F  layer  alti¬ 
tudes.  Few  electrons  had  sufficient  energy  to  create  E-region 
ionization  (Figure  6c).  Recall,  too,  that  in  the  evening 
sector,  precipitating  ions  carried  the  downward  FAC 
(Figure  6d).  Galand  and  Richmond  [2001]  and  Galand  et 
al.  [2001]  showed  that  precipitating  protons  significantly 
affect  ionospheric  conductance.  Determining  quantitative 
relationships  between  EP  and  precipitating  ions/electron 
fluxes  is  beyond  the  scope  of  this  report. 

[61]  Uncritical  applications  of  the  Robinson  formulas, 
irrespective  of  location  on  the  globe,  can  also  lead  to 
significant  errors.  Assume  that  at  6BZ  extrema,  the  SSJ4 
sensors  on  F13  and  FI 5  measured  the  same  values  of  $E 
and  EAVe-  The  Robinson  formulas  predict  the  same  values 
of  EP  and  EH  beneath  the  two  satellites.  This  cannot  be 
correct.  Pedersen  conductivity  is  given  by  crP  «  (ne/B)  • 
(■Ojn/wj),  where  n,  uin,  and  W;  represent  the  plasma  density,  the 
ion-neutral  collision  frequency,  and  the  ion  gyrofrequency, 
respectively.  To  a  first  approximation,  electrons  precipitat¬ 
ing  at  the  locations  of  F13  and  FI 5  would  produce  the  same 
plasma  density  profiles.  Assuming  similar  thermospheric 
profiles,  the  uin  height  profiles  should  also  be  similar. 
However,  since  o>  is  proportional  to  B~2  in  the  Pedersen 
current  layer,  values  of  EP  =  f  aP  ds  should  be  quite 
different.  The  magnetic  field  measured  by  F13  was  1.58 
times  that  observed  by  FI 5.  Thus  for  the  same  particle 
inputs,  EP  at  the  location  of  FI  3  would  be  ~40%  of  that  at 
the  location  of  FI  5.  At  the  interface  between  the  two  FACs, 
we  estimated  that  at  the  location  of  F13  EP  was  46%  of  the 
value  at  F15.  On  this  empirical  basis,  we  join  Kamide  et  al. 
[1981]  and  Raeder  et  al.  [2001]  in  cautioning  that  integra¬ 
tion  of  space  measurements  into  AMIE  simulations  of 
superstorms  must  be  exercised  with  great  care  less  they 
produce  misleading  results. 

5.2.3.  Transmission  Line  Effects 

[62]  From  an  electrodynamic  point  of  view,  obliquely 
propagating  AlfVen  waves  carry  the  FACs  [Kan  and  Sun, 


1985],  Quasi-steady  state  currents  represent  standing  wave 
structures  resulting  from  the  superposition  of  incident  (i) 
and  reflected  (r)  waves.  Electric  and  magnetic  field  sensors 
on  polar-orbiting  satellites  detect  superposed  perturbation 
fields.  In  the  infinite  current  sheet  approximation,  the 
fields  are  EY  =  E'Y  +  EY,  and  6BZ  =  6B^  +  6Brz. 
The  incident  and  reflected  electric  fields  are  related 
through  the  relationship  EY  =  R  E'Y,  where  the  reflection 
coefficient  R  =  (EA  —  EP)/(EA  +  EP).  The  term  EA  = 
l/(poVAR)  is  called  the  AlfVen  conductance;  VAR  is  the 
Alfven  speed  in  the  reflection  layer.  Observed  above  the 
reflection  layer,  the  incident  and  reflected  waves  are  seen 
as  propagating  at  the  local  Alfven  speed  ±VAS,  i.e., 
E'y/&Bz  =  VAS  =  —  Ey/6BZ.  Thus  &BZ  =  -R  6Brz. 

[63]  The  ratio  of  magnetic  to  electric  field  variations  at 
the  satellite  altitude  is  given  by 

Wl  =  Wz  +  bBrz  =  1  1  -R=  1  E  P 

Er  \  measured  =  E‘Y+ErY  ( Vas)1+R  (Vas)Xar'  {  } 

The  IGRF  model  provides  estimates  of  B  at  any  altitude  in 
the  ionosphere.  Local  plasma  densities  are  measured  by 
SSIES.  Across  regions  of  interest,  0+  ions  were  the 
dominant  species.  With  this  information,  we  can  calculate 
local  values  of  the  Alfven  speed,  VAS.  Knowing  VAR 
requires  unavailable  information  about  the  density  and 
mean  mass  of  ions  in  the  reflection  layer.  However, 
equation  (16)  allows  us  to  use  measured  values  of  VAS, 
6BZ,  and  Ey  to  calculate  the  ratio  EP/EAR  and  determine  the 
reflection  coefficients  R. 

[64]  Table  3  combines  6BZ  and  Ey  observed  by  F13  and 
FI 5  in  upward  and  downward  FAC  sheets  (Table  2)  with 
measured  values  of  VAS  to  calculate  EP/EA  and  R. 
Columns  5  and  6  of  Table  3  show  the  calculated  intensities 
of  electric  and  magnetic  fields  for  incident  Alfven  waves 
Ey  =  Ey/(1  +  R)  and  6BZ  =  6BZ/(1  —  R).  Except  for  the 
downward  FAC  crossed  by  FI 3,  reflection  coefficients  are 
near  0.9.  In  all  cases  the  Ey  values  are  much  larger  than  are 
typically  observed  at  840  km.  Note,  however,  that  in  all  cases 
the  ratio  AEY/A6BZ  agrees  with  the  local  Alfven  speed 
derived  from  plasma  densities  and  total  magnetic  field 
measurements.  Poynting  fluxes  carried  by  the  incident  waves 
are  of  order  0.1  W/m2  =  1 00  erg/(cm2  s).  Although  50  to  80% 
of  this  energy  is  reflected,  the  Joule  heat  absorbed  by  the 
ionosphere  and  thermosphere  greatly  exceeds  the  energy 
deposited  by  precipitating  electrons  and  ions. 

[65]  Our  analysis  of  DMSP  electric  and  magnetic  field 
measurements  provides  information  about  the  ratio  EP/EA 
but  nothing  about  the  value  of  either  quantity.  Data  presented 
in  the  third  column  of  Table  3  indicate  that  EP/EA  ratios  are 
relatively  large.  At  ionospheric  altitudes  Alfven  impedances 
are  typically  of  the  order  of  1  mho.  Thus  EP  values  are  easy  to 
reconcile  with  estimates  in  Table  2.  On  the  other  hand,  if  we 
assume  that  EP  has  values  determined  using  the  formulas  of 
Robinson  et  al.  [1987],  then  proportionately  lower  values  of 
Ear  are  required.  This  in  turn  requires  that  VAR  >  VAS,  a 
conclusion  that  is  difficult  to  reconcile  with  plasma  density 
profiles  that  typically  decrease  exponentially  with  altitude  in 
the  topside  ionosphere. 

[66]  The  fact  that  SSIES  and  SSM  detect  superposed 
downward  and  upward  propagating  Alfven  waves  has 
practical  consequences.  Measurements  of  6BZ  and  EY  can 
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be  used  to  estimate  the  Poynting  fluxes  associated  with  the 
current-carrying  waves.  Equation  (10)  indicates  that  the  net 
Poynting  flux  is  a  direct  measure  of  Joule  heat  dissipated  in 
the  ionosphere  magnetically  conjugate  to  the  spacecraft. 
Attention  is  redirected  to  the  large  FAC  structure  centered 
near  203 1  UT  in  panel  4  of  Figure  5c.  We  have  calculated 
Poynting  fluxes  between  55.5°  and  61.7°  MLat.  The  max¬ 
imum  net  Poynting  flux  was  ~0.11  W/m2.  The  integrated 
electromagnetic  energy  flowing  into  the  ionosphere  across 
the  structure  was  ~4.2  x  104  W/m  (cf.  equation  (11)). 
Measurements  from  multiple  DMSP  satellites  indicate  that 
the  large  FAC  structures  extend  for  at  least  an  hour  in  local 
time.  Over  a  15°  spread  in  longitude  the  total  power  input  to 
the  ionosphere  exceeds  0.04  TW.  DMSP  data  also  indicate 
that  the  lifetime  of  these  structures  is  ~20  min. 

[67]  If  we  assume  that  the  FAC  systems  detected  by 
DMSP  satellites  occurred  in  both  hemispheres  and  extended 
for  3  hours  in  local  time  in  both  the  evening  and  morning 
sectors,  their  total  Joule  heat  input  to  the  ionosphere  in 
20  min  would  be  ~  144  TJ.  This  represents  ~4%  of  the 
total  kinetic  energy  of  the  ring  current  particles  estimated 
using  the  Dessler-Parker-Sckopke  relation  with  Dst  = 
—250  nT  [Carovillano  and  Maguire,  1968].  We  suggest 
that  FACs  can  act  as  significant  coolants  for  ring  current 
populations  and,  along  with  intense  fluxes  of  precipitating 
ions  [Borovsky  et  al.,  1998],  partly  explain  differences 
between  predicted  and  observed  ring-current  ion  spectra 
during  the  June  1991  magnetic  storm  [Garner  et  al.,  2004]. 

[68]  Finally,  our  analyses  relating  FACs  to  Alfven  waves 
and  to  magnetospheric  pressure  gradients  are  not  necessarily 
self-consistent.  The  Vasyliunas/RCM  model  describes  quasi¬ 
steady  state  developments  of  magnetospheric  electrodynam¬ 
ics  that  are  averaged  over  several  Alfven  bounce  periods. 
Song  and  Lysak  [2000]  have  pointed  out  that  current- 
carrying  Alfven  waves  can  be  generated  via  the  mode 
conversion  of  compressive  magnetohydrodynamic  (MHD) 
waves.  Analyzing  the  possible  role  of  MHD  mode  conver¬ 
sion  is  well  beyond  the  scope  of  this  observational  report. 
However,  we  note  that  the  compressive  modes  are  intimately 
tied  to  variations  of  magnetic  field  intensities  and  plasma 
pressures.  Our  superstorm  data  suggest  that  the  growth  and 
decay  of  large  FAC  structures  in  the  ionosphere  reflect  the 
dynamics  of  magnetospheric  plasma  pressure  distributions. 
Uncovering  the  physical  details  of  mechanisms  responsible 
for  large  amplitude,  ultra  low  frequency  Alfven  waves 
associated  with  FACs  represents  a  fruitful  avenue  in  space 
research  to  expand  understanding  of  M-I  coupling. 

5.3.  Absence  of  Ground  Magnetic  Responses 

[69]  In  absolute  terms  ground  magnetometers  continued 
to  observe  significant  perturbations  during  the  interval  after 
2000  UT.  SYM-H  and  AE  indices  were  near  —200  and 
500  nT,  respectively  (Figure  2).  Ring  current  effects  are 
clearly  visible  in  magnetic  records  from  stations  at  low  to 
middle  magnetic  latitudes.  Relatively  high  levels  of  activity 
persisted  at  the  latitudes  of  auroral  electrojet  stations,  well 
poleward  of  the  transient  but  intense  FAC  sheets  reported 
here.  At  first  glance  the  fact  that  FACs  associated  with 
|6Bz|  >  1200  nT  produce  ground  magnetic  perturbations 
<100  nT  appears  surprising.  We  note,  however,  that  near 
the  extrema  of  6BZ  and  EY  values  of  £H  estimated  from 
fluxes  of  precipitating  electrons  are  very  low  (Figure  9). 


Fukushima  [1976]  predicted  that  FACs,  closed  only  by 
Pedersen  currents,  produce  no  magnetic  perturbations  at 
the  ground.  Perturbations  detected  on  the  ground  mostly 
come  from  Hall  currents.  With  low  values  of  £H,  Hall 
currents  are  weak.  We  contend  that  the  combined  DMSP 
and  ground  measurements  constitute  an  empirical  validation 
of  Fukushima’s  prediction. 

5.4.  Precipitation  Fluxes 

[70]  As  mentioned  above,  the  dominant  species  observed 
by  DMSP  F12  in  the  interval  from  2021:15  to  2022:30  UT 
appear  to  be  precipitating  ions.  This  can  be  seen  from  the 
number  flux  plot  in  the  panel  2  of  Figure  5b.  Figure  6 
explicitly  shows  the  difference  between  electron  and  ion 
number  fluxes.  In  the  upward  FAC  region,  there  is  a  net  flux 
of  precipitating  electrons.  However,  in  the  downward  FAC 
region  the  precipitating  ion  flux  dominates.  Ion  and  electron 
spectra  on  the  right  side  of  Figure  6  illustrate  the  difference 
in  the  particle  populations  with  access  to  the  ionosphere  in 
the  bipolar  current  sheets. 

[71]  The  SSJ4  measurements  raise  fundamental  questions 
for  modelers:  What  dynamics  are  responsible  for  episodes  in 
which  normally  trapped  ions,  over  a  broad  range  of  energies, 
are  suddenly  dumped  into  the  ionosphere?  Examining  ion 
spectra  at  magnetic  latitudes  between  B2e  and  B2i  measured 
by  F12  and  F14  before  and  after  the  FAC  event  shown  in 
Figure  6  leads  us  to  believe  that  the  ion  deposition  of  Figure 
6d  represents  a  temporal  effect.  A  broad  spectrum  of  ions  that 
were  previously  trapped  on  flux  tubes  mapping  to  MLat 
»55°  precipitate.  The  energy  range  of  affected  ions  seems  to 
rule  out  pitch-angle  scattering  via  gyroresonant  wave-particle 
interactions  [Kennel  and  Petschek,  1966],  In  the  magneto¬ 
sphere  trapped  ions  follow  guiding-center  trajectories  along 
which  the  first  (mv2L/B)  and  second  (J  m;V||  ds)  adiabatic 
invariants  are  constants  of  the  motion.  If  the  length  of 
affected  magnetic  flux  tubes  decreases  during  intervals 
significantly  longer  than  an  ion  bounce  period,  in  conser¬ 
ving  the  second  invariant  a  particle’s  Vy  increases  and  its 
pitch  angle  decreases.  This  lowers  mirror  altitudes,  allowing 
some  previously  trapped  ions  to  become  lost  in  the 
ionosphere.  Why  would  such  a  field-line  shortening  occur? 
One  possible  source  is  a  compression  of  the  magnetic  field 
due  to  an  injection  of  plasma  from  the  magnetotail.  This 
is  not  the  place  to  answer  these  questions  raised  by  the 
DMSP  measurements.  However,  from  either  quasi  steady- 
state  [ Vasyliunas ,  1970]  or  dynamical  [Song  and  Lysak, 
2000]  points  of  view,  the  answers  involve  interactions 
between  plasma-pressure  and  magnetic-field  distributions  in 
the  highly  stressed  inner  magnetosphere  that  are  still  not 
understood. 
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MLT  2042  20.7.1  20. '>5  21.13  19  5  19  X  20  0  20.13  07.24  <16  X4  06.68  06.55 

Figure  5.  Particle  and  field  measurements  by  (a)  DMSP  FI 5  and  (b)  FI 2  in  the  evening  sector  after 
2015  UT,  and  (c)  FI 3  in  the  morning  sector  after  2024  UT  on  6  April  2000.  From  top  to  bottom  the  panels 
for  each  satellite  show:  (1)  directional  energy  flux  (E-Flux)  in  eV  cm-2  s-1  sr-1  of  down-coming 
electrons  (upper  line)  and  ions  (lower  line);  (2)  directional  number  fluxes  (N-Flux)  of  down-coming 
electrons  (upper  line)  and  ions  (bottom  line)  with  energies  between  30  eV  and  30  keV,  in  cm~2  s-1  sr~  , 
(3)  average  energy  of  precipitating  ions  (upper  line)  and  electrons  (lower  line),  defined  as  EAVe  = 
E-Flux/N-Flux;  (4)  the  cross-track  component  of  the  magnetic  field  8BZ  in  nanotesla;  and  (5)  the  in-track 
component  of  the  electric  field  EY  calculated  from  the  observed  ion  drift  and  IGRF  magnetic  field  values 
[Rich  and  Hairston,  1994].  Vertical  dashed  lines  indicate  times  when  the  DMSP  satellites  crossed  Bl,  B2e, 
B2i,  and  B5  [ Newell  et  al.,  1996],  Data  are  plotted  as  functions  of  UT,  MLat,  and  MLT. 
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Figure  6.  Particle  and  field  characteristics  of  the  FAC  sheets  crossed  by  DMSP  F12  between  2020  and 
2023  UT  on  6  April  2000.  Panels  show  (a)  variations  of  8BZ  and  (b)  differences  between  electron  and  ion 
number  fluxes  as  a  function  of  UT,  as  well  as  5-s  averaged  ion  and  electron  spectra  measured  by  SSJ4 
within  the  (c)  upward  and  (d)  downward  FAC  sheets. 
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14.  ABSTRACT 

[t]  During  the  main  phase  ot  the  6  April  2000  magnetic  storm  with  Dst  =  —300  nT  four 
Defense  Meteorological  Satellite  Program  (DMSP)  satellites  encountered  intense  sheets  of 
field-aligned  currents  (FACs).  Their  magnetic  perturbations  were  >1300  nT 
corresponding  to  integrated  currents  |J„|  >  1  A/m.  The  FACs  appeared  in  both  the  evening 
and  dawn  magnetic  local  time  sectors.  They  had  relatively  fast  rise  times  (~5  min),  lasted 
for  ~0.5  hours,  and  were  associated  with  widespread  reconfigurations  of  plasma  in’ the 
near-Earth  magnetotail.  The  largest  magnetic  and  related  electric  field  perturbations 
occurred  at  magnetic  latitudes  <60°.  Magnetometer  measurements  from  DMSP  satellites 
show  repeated  episodes  of  similarly  large  FACs  late  in  the  main  phase  of  this  and 
'  other  superstorms.  Poynting  flux  calculations  indicate  that  a  few  percent  of  the  total  ring 

current  energy  is  dissipated  as  Joule  heat  in  the  midlatitude  ionosphere  during  each  of 
these  events.  . 
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